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We have introduced recently a model for the evolution 
of a sexual transmitted disease where the social behavior 
is incorporated by means of what we call the promis- 
cuity variableQ, i.e. the daily probability of an agent 
to dismiss its permanent mate and go out to look for a 
sexual intercourse. The dynamics is of the SIR type: a 
susceptible (S) agent can become infected (I) after a sex- 
ual intercourse with an infectious one, with probability 
f3 (the per sexual contact probability of get infected by 
an infectious agent). The infected individuals remain in- 
fectious for a period r, when they are removed (R) by 
death. The system will model the dynamic network of 
social contacts, where each element of the network is con- 
sidered as a sexually active subject. The social behavior 
element is included via the promiscuity and the fraction 
of singles. In principle, individuals are grouped in couples 
and additionally a controlled proportion of singles, (ps), 
is considered. The variable pi, called the "promiscuity", 
is what determines an individual's tendency to dismiss 
his/her mate and go out -or just to go out in the case of 
the singles- to look for an occasional intercourse; more 
precisely pi is the probability of trying to meet an oc- 
casional partner on each intercourse opportunity or time 
step of the simulation (which, for the sake of simplicity, 
we take as a day). The pi's are randomly taken from a 
chosen distribution and assigned to each individual; we 
start with a semi-Gaussian distribution of width equal to 
p, but we will compare later others distributions. Those 
who happened to go out chose a partner at random and if 
the latter happened to go out too, the occasional couple 
is made. Therefore there is no social structure a priori. 
The web of contacts is constructed dynamically during 
the simulation depending on the percentage of singles 
and the population promiscuity. Most of the results are 
for the homosexual case, but we considered the hetero- 
sexual case too. The details of the simulation are given 
in Q , and the main result of this model can be expressed 
in the following relation, relating the four parameters of 
the model at the epidemic transition: 



respectively for 



plf3T = 0.36,0.48,0.65 



P. = 0,0.5,1 



where p is a per day probability, /? is a per intercourse 
probability, and t is in days. This relation gives us the 
following values for the mean promiscuity at the critical 
point: 0.025, 0.029, and 0.034, respectively for the three 
values of ps, and for/? = 1 and r = lyr. 

Some authors Q y, Q have pointed to the connection 
between epidemics and percolation. In our case there 
is no network a priori, no topology with fixed links be- 
tween nodes (agents); instead, the sexual network is con- 
tinuously built and destroyed as time goes by, depending 
upon the individual promiscuity, that are controlled by 



the promiscuity distribution, and the fraction of singles. 
Therefore the threshold can be understood as a dynamic 
percolation. This is supported not only by the sharp 
transition of the asymptotic number of death at the crit- 
ical promiscuity value, but by the divergence of the time 
needed to get that asymptotic number of death (remem- 
ber that there is no other cause of death in the present 
model) at the same promiscuity value, as can be seen 
in Fig. All of this has a strong resemblance with a 
percolation transition. 
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FIG. 1: Relaxation time as a function of average promiscuity 
for ps = 0, for /3 = 1 and r = lyr. 

In the classical implementation of the SIR model 
(see Q), the dynamics of the susceptible, infectious, and 
removed classes is governed by the following ordinary 
nonlinear equations: 



ds 
'dt 



di 
dt 



ris 



dr 
'dt 



ai 



where r is the infectious rate and a is the removal rate 
per unit time; r can be expressed as the product of (3 and 
c, the mean number of sexual contact per unit time, and 
a is the inverse of the infectious time r. The onset of 
the epidemic happens when the basic reproductive num- 
ber i?o = fire is equal to 1. Therefore it is interesting 
to compare the threshold and the number of asymptotic 
death of the standard fully mixed model with the present 
results. In the former case we have to resort on numerical 
integration to solve the equations above. We do this up 
to a time were no new susceptible get infected (in prac- 
tice when dS < l/N) and we plot the final number of 
removed agents against the average number of individual 
sexual contact per unit time, assuming that 13—1 and 
T = lyr. The results are plotted in Fig 12 together with 
the prediction of our model. In order to do the compar- 
ison we have obtained from the simulations the relation 
between the average promiscuity and the average num- 
ber of contacts per unit time (the latter depends almost 
quadratically on the former one). 
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FIG. 2: Asymptotic percentage of removed agents as a func- 
tion of the average number of contacts per year, for three 
different values of ps Comparison with the numerical solution 
of the fully mixed SIR model (/3 = 1 day~^ , t — 1 yr). Sim- 
ulations are made with 100000 subjects and 100 realizations. 



Comparing them we see that the threshold predicted 
by the numeric solution of our model is considerably 
lower, possibly due to the effect of fluctuation that the 
uniformly mixed SIR model does not take into account. 
Moreover, the divergence between the latter and simu- 
lations increases when couples got into the scene. This 
more sophisticated version speed up the epidemic so 
we conclude that the epidemic dynamics, in general, can 
not be reduced to take the average number of contacts 
as the only relevant parameter; as for our results, ps for 
example, is a relevant parameter too, and all of this re- 
inforce the idea of using simulation as the appropriated 
technique to study epidemics. 
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FIG. 3: Asymptotic percentage of removed agents as a func- 
tion of the mean promiscuity of the population (p) for sev- 
eral promiscuity distributions. Parameters are: r = lyr and 
P = lday~^ . Numerical results obtained with 100000 sub- 
jects, 100 realization, and for ps — 1 case (all singles). 



The results summarized in Fig. and Fig. 12 are for a 
semi-Gaussian distribution of promiscuity, so the natural 
question is "How sensitive to the choice of distribution 
are the results presented here?" In Fig. we plot the 
results for a delta {pi = constant) and a uniform distri- 
bution together with the semi-Gaussian distribution, all 
of them normalized to the corresponding mean promis- 
cuity {p). We can see in Fig. El that the uniform and 
the semi-Gaussian distributions present almost the same 
threshold, while the delta distribution has its threshold 
slightly shifted to a larger promiscuity, due to the absence 
of heterogeneity. In the other hand we have previously 
shown that a power law distribution -more appropri- 
ated to the real word [^-, produce a contrary effect: a 
threshold shifted to a smaller promiscuity value. Never- 
theless even in this latter case it is a weak effect, so we 
conclude that irrespective of the distribution of promiscu- 
ity, that in turn governs the distribution of contacts, the 
mean value is the most relevant measurements to mark 
the onset of the epidemic, while other moments (espe- 
cially the standard deviation) seem not to be important, 
in contradiction with recent discussions 0, 

A flnal word about the solution of the fully mixed clas- 
sical implementation of the SIR model: generally speak- 
ing it is said that the limitation of this approach came 
from the fully mixed assumption and from assigning the 
same value of infectivity, and the same number of per 
year contacts to all individuals. However, apart from 
that limitations, there is one that comes flrst, it cames 
from the method of solution that assumes the population 
as continuous quantities without fluctuations. Let's see: 
if we take the delta distribution, i.e. the same promis- 
cuity for all agents, we have the model in the condition 
described above: each agent takes a potential partner at 
random from the whole population, so it is indeed fully 
mixed, and the parameters p, /?, and r are the same for 
everybody. However if we compare the curve of the clas- 
sical SIR model solution of Fig. |21 with the results for 
the delta distribution in Fig. the differences are ob- 
vious. Not only the threshold is overestimated by the 
classical approximation, the simulation predicts a sharp 
transition separating the no epidemic region from the epi- 
demic one with no survivors at all. We emphasize that, 
apart from the limitation of the model, we are comparing 
here the exact solution of it (simulations) with an con- 
tinuous -approximated- solution of the same model (the 
ODE equations) 

So far we have assumed an heterogeneous population 
in terms of the promiscuity, but homogeneous on all of 
other respects, more appropriated for a gay community 
for example. What happens if we consider a heterosexual 
population? Taking the all married case version, but now 
50% males and 50% females, we have the same output 
as for the homosexual case, as can be observed in Fig.^ 
The point is that even when we have two subgroups in the 
population there is no asymmetry (the same promiscuity 
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FIG. 4: Asymptotic percentage of removed agents as a func- 
tion of mean population promiscuity (p). Numerical result 
obtained with 100000 subject and 100 realization for ps = 
case. Heterosexual and homosexual cases are displayed with 
symbols for each realization; the line is an average over them. 

distribution, infectivity and,r for all agents, whether they 
are males or females) so the results is no surprising. If we 
broke the symmetry we would expect a different behavior 
between the hetero and homo cases as we plan to study 
in the near future. 

Summarizing we have presented some results for the 
spread of sexual transmitted diseases, based on a model 
that takes into account the social behavior. The model 
predicts a threshold condition quite different that what 
the classical SIR models predict. The onset of the epi- 
demic is associated in the present model with a perco- 
lation threshold of a dynamical kind, because there is 
no static network of contact. The present model gives a 



direct connection between disease and social parameters 
and the outbreak threshold, which is quantitative equiv- 
alent to a basic reproductive number Rq = 0.64, for the 
singles case (0.4 for the all married case), much lower 
than that of the fully mixed model prediction (i?o = !)• 
On the other hand we have shown that the onset of epi- 
demic depends mainly on the average number of contact, 
but seems to be independent of the exact form of the dis- 
tribution, in opposition to recent analysis that include 
the second moment in the basic reproductive number of 
heterogeneous populations. Finally the threshold is sen- 
sitive to the fraction of singles, and this can have impor- 
tant consequences in modeling real data. 
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